Local atomic structure and discommensurations in the charge density wave of CeTe 3 
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The local structure of CeTe3 in the incommensurate charge density wave (IC-CDW) state has 
been obtained using atomic pair distribution function (PDF) analysis of x-ray diffraction data. Local 
atomic distortions in the Te-nets due to the CDW are larger than observed crystallographically, 
resulting in distinct short and long Te-Te bonds. Observation of different distortion amplitudes 
in the local and average structures are explained by the discommensurated nature of the CDW 
since the PDF is sensitive to the local displacements within the commensurate regions whereas the 
crystallographic result averages over many discommensurated domains. The result is supported by 
STM data. This is the first quantitative local structural study within the commensurate domains 
in an IC-CDW system. 

PACS numbers: 71.45.Lr, 61.44.Fw, 61.10.Nz 



Incommensurate charge density waves (IC-CDWs) are 
a fundamental property of low-dimensional metals and 
also underly the novel properties of correlated electron 
oxides such as cuprates in the pseudo-gap state d El 
and manganites at high doping 0]. Knowing the na- 
ture of local atomic displacements (Peierls distortions) 
in the IC-CDWs is crucial to understand such factors 
as electron- lattice coupling yet this information is 
difficult to obtain quantitatively. Here we solve this 
problem by taking the novel approach of using a local 
structural method, the atomic pair distribution function 
(PDF) technique [jj, to determine the local atomic dis- 
placements with high precision in the system CeTe3. IC- 
CDWs, and the underlying atomic displacements, can 
be uniform incommensurate modulations or locally com- 
mensurate waves separated by narrow domain walls, 
known as discommensurations [8j , where the phase of the 
wave changes rapidly. Here we show that the IC-CDW 
in CcTe3 is discommensurated and obtain for the first 
time the quantitative local atomic displacements within 
the commensurate domains. 

In the case of incommensurate CDWs, superlatticc 
peaks observed crystallographically yield the average dis- 
torted structure. Except in the cases where the domains 
are periodically arranged, giving rise to satellite peaks , 
it is not possible to determine whether the underlying 
CDW is truly incommensurate or forms a discommen- 
surated structure with commensurate regions separated 
by domain walls Q. A number of techniques have been 
successful at differentiating between the truly incommen- 
surate and discommensurated cases. The earliest ver- 
ification of a discommensurated phase came from pho- 
toemission spectroscopy evidence that the Ta 4/ states 
in lT-TaS2 had the same splitting in the commensu- 
rate and nearly-commensurate states [T^j. Photoemis- 
sion is a local probe and found distinct Ta environ- 
ments rather than a broad continuum expected from 
a purely incommensurate state. Similarly, another lo- 



cal probe, nuclear magnetic resonance (NMR), found 
distinct Knight-shifts for three Se sites in the incom- 
mensurate state of 2H-TaSe2, similar to the commen- 
surate phase 0, 0| . High resolution atomic imag- 
ing methods have also contributed to this debate. The 
strain fields due to the domain walls were observed in 
dark field transmission electron microscopy (TEM) mea- 
surements [lflj . Interestingly, atomic resolution images 
in real-space have difficulty in resolving discommensu- 
rated domains [Til Hoi flil Il7j . However, Fourier analy- 
sis of scanning tunneling microscopy (STM) images can 
be a reliable measure, as discussed in detail by Thom- 
son et al. [Isj . 

As in the case of the NMR and photoemission studies, 
the PDF approach described here makes use of the fact 
that the local structure deviates from the average in the 
discommensurated case. By comparing atomic displace- 
ments determined from the PDF with those determined 
crystallographically we establish the presence of commen- 
surate domains, but crucially, also obtain quantitatively 
the atomic structure within these domains. This novel 
approach is here applied to the incommensurate phase of 
CeTe 3 . 

In its undistorted form, CeTe3 takes the NdTe3 struc- 
ture type with space group Cmcm . It forms a layered 
structure with ionic [Ce2 3+ Te2 2_ ] 2+ layers sandwiched 
between two Te~ layers. These sandwich layers stack 
together with weak van der Waals forces to form the 3- 
dimensional structure. Te ions in the Te - layers form 
a square-net with 3.1 A Te-Te bonds. The structure is 
shown in Fig. H{a). The electronic bands crossing the 
Fermi level are Te p-bands from the 2D square nets [^ 
and the CDW forms in these metallic layers. In the CDW 
state an incommensurate superlattice is observed |2lj . 
with a wavevector characteristic of a stron g Fe rmi-surface 
nesting vector in the electronic structure [22I I23L I24I E^f . 
This is a surprisingly stable and simple single-g IC-CDW 
state in an easily cleavable 2D square net making the 
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FIG. 1: (a) The crystal structure of CeTe3 with the square 
Te net that supports the CDW highlighted. The reduced unit 
cell on the Te net is indicated by the red dashed box (- -). (b) 
A representative STM image from the square Te net showing 
the CDW. On the expanded image, the network of Te bonds 
is superimposed, (c) The Fourier transform of the STM data. 
To achieve a high signal-to-noise ratio, the transform repre- 
sents the average of 24 images (each image was 27 x 27 nm). 
The square Te net gives rise to four distinct peaks (L), with 
peaks related to the CDW oriented at 45°, as indicated by 
the arrow. The fundamental CDW peak (corresponding to a 
wavelength of w 15 A) and the A/2 harmonic are labeled 1 
and 3, respectively. Peaks 2 and 4 are in close proximity to 3, 
implying a characteristic discommensuration length of 38 A, 
as described in the text. Peak 5 corresponds to the diagonal 
of the Te net. This component of the Te net may be enhanced 
due to the underlying crystal structure; the CDW-lattice in- 
teraction may also enhance this peak. 



RETC3 (RE=Rare Earth) systems ideal for studying the 
IC-CDW state [13 • The atomic distortions giving rise 
to the superlattice have been solved crystallographically 
from single crystal x-ray diffraction data |2lj . The incom- 
mensurate wavelength of the distortion is close to 25a/7, 
where a is the lattice parameter of the undistorted phase. 
The distorted structure is in the Ama2 spacegroup [2lj . 
From the crystallography alone it is not possible to de- 
termine whether this distorted structure is truly incom- 
mensurate or whether discommcnsurations form between 
short-range commensurate domains. 
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FIG. 2: (a) The PDF of CeTe 3 at room temperature. In 
(b)-(d) the first peak of the experimental PDF of CeTe3 (o) 
is plotted with the calculated PDF (— ) from various models: 
(b) the undistorted crystal structure model (Cmcm), (c) the 
distorted crystallographic model and (d) the local structural 
model determined by the PDF refinement over the range 2.5 < 
r < 6.37 A. The difference between the experimental and 
calculated PDFs (— ) is plotted below the data in each panel. 
The shoulder due to the Peierls distortions in the Te-nets is 
indicated by an arrow. 



The X-ray PDF experiment was conducted on a fine 
powder of CeTe3 prepared as described in Ref. [ilj . 
CeTe3 powder was loosely packed in a flat plate with 
thickness of 1.0 mm sealed with kapton tape. Care must 
be taken when grinding this material or turbostratic dis- 
order is introduced, significantly modifying the stacking 
of the layers. Diffraction data were collected at 300 K 
using the rapid acquisition pair distribution function 
(RA-PDF) technique (H- Standard corrections 0, El 
were made using the program PDFgetX2 [2^ to ob- 
tain the properly normalized total scattering function, 
S{Q), which was truncated at Q ma x of 25 A^ 1 be- 
fore Fourier transforming to obtain the PDF, G(r) = 
f Jo°° Q[S(Q) — 1] sin(Qr) dQ. Structural models are fit 
to the data using the program PDFFIT |2s| . 

The PDF of CcTc3, measured at room temperature, 
is shown in Fig. |2{a). The PDF gives the probability 
of finding an atom at a distance-r away from another 
atom. The nearest neighbor peak around 3.1 A comes 
from the Te-Te bond in the nets and the Ce-Te bond in 
the intergrowth layers. This is shown on an expanded 
scale in Figs.|2^b)-(d). A shoulder is evident on the low- 
r side of the peak. This feature is robust; it is much 
larger than the statistical and systematic errors and is 
reproduced in measurements of isostructural compounds 
NdTe 3 and PrTc 3 . Fig.^b) shows the fit to this peak of 
the undistorted crystal structure model (Cmcm), where 



3 



only symmetry allowed atomic positions and isotropic 
thermal factors were refined. The result clearly does 
not explain this shoulder which originates from short Te- 
Te bonds in the Te-net. Surprisingly however, the PDF 
calculated from the distorted structure determined crys- 
tallographically [2lJ also does not explain this shoulder 
well. In this case the atoms were fixed at the crystal- 
lographically determined positions and isotropic thermal 
factors were refined. This resulted in a better fit to the 
first peak (Fig. EJc)); however, the fit is not ideal and 
required a large value of Ui SO for the Te atoms in the 
nets (U iso = 0.0152(2) A 2 ). The value was two times 
larger than Ui SO of the Ce and Te atoms in the ionic 
[Ce 2 3+ Te 2 2 -] 2+ layers (0.0077(2) A 2 and 0.0080(2) A 2 
for Uiso of Ce and Te atoms, respectively). 

The large fluctuation in the difference curve in 
Fig. |2fb) arises because the real distribution of Te-Te 
bond-lengths in the data is broader than in the undis- 
torted model. This fluctuation in the difference curve is 
smaller in Fig.EJc) because the distortions of the Te net 
in the Ama2 crystallographic model result in a broader 
Te-Te bond-length distribution. However, clearly the 
distorted-crystallographic model still has a Te-Te bond- 
length distribution that is narrower than in the data. 
We therefore refined the Te-net distortions directly in 
the PDF by allowing the atomic positions in the model 
to vary. The model was constrained to have the Amal 
symmetry and the same unit cell was used as in the dis- 
torted crystallographic model. The refinement result for 
2.5 < r < 6.37 A is shown in Fig. E(d). As well as giv- 
ing a significantly better fit to the low-r region of the 
PDF, this refinement resulted in much smaller and more 
physical thermal factors on the planar Te ions. 

The model of the local structure refined from the PDF 
gives a broader range of Te-Te bond lengths (from 2.83 A 
to 3.36 A) than the crystallographic distorted model 
(from 2.94 A to 3.26 A). It is also interesting to see 
the shape of the bond-length distributions for the Te- 
Te bonds in the Te-nets from these two models. These 
are shown in Fig. El T ne blue solid line shows the bond- 
length distribution refined from the PDF and the red line 
is the distribution from the crystallographic model [2l| . 
For direct comparison the distributions are plotted using 
the same thermal broadening of 0.007 A 2 . The distorted- 
crystallographic model has broad, but symmetric and 
Gaussian bond-length distribution coming from the con- 
tinuous distribution of Te-Te bond lengths in the aver- 
age IC-CDW. On the other hand, the local structure re- 
finement (r max = 6.37 A) yields a bond-length distribu- 
tion that is clearly bimodal and is separated into distinct 
"short" and "long" Te-Te distances. This is emphasized 
in Fig. |^b) where we show a fit of two, well separated, 
Gaussian curves to the PDF-refined bond-length distri- 
bution. This behavior is characteristic of oligomerization 
with Te forming bonded and non-bonded interactions 
with its neighbors in the net [2^ that would be expected 
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FIG. 3: Bond- length distributions in the Te-nets refined from 
the PDF over various r-ranges (— ) (a) r maa; =6.37 A (b) 
r ma x=6.37 A (c) r roMI =14.5 A (d) r max =27.1 A. The r min 
value was fixed to 2.5 A for all the cases. The bond distribu- 
tion from the distorted crystallographic model (— ) is plotted 
in (a), (c), and (d) for comparison. In (b) the bond distribu- 
tion of the local structural model (r max —6.37 A) is fit with 
two Gaussians. The fit is shown as a red line (— ) and the two 
Gaussian sub-components in green (— ). 



in a commensurate structure. Since we know that the 
modulation is incommensurate on average, this is strong 
evidence that the structure consists of commensurate do- 
mains separated by discommensurations. As r max in the 
PDF-refinements is increased the PDF refined distribu- 
tion crosses over towards the crystallographic result and 
by r max = 27.1 A, resembles it rather closely (Fig. EJd)) • 
We have applied STM on the exposed Te net of a 
cleaved single crystal of CeTe3, grown according to the 
method described in Ref. [2l|. Measurements were done 
at 300 K in the constant current mode of the STM. Data 
were acquired with a bias voltage of 100 mV and with a 
tunneling current of 0.6 nA. Fig.^b) shows a representa- 
tive atomic resolution image with the CDW modulation 
clearly visible oriented at 45° to the net. To investigate 
the images for discommensurations, we examine the cor- 
responding two-dimensional Fourier transform, shown in 
Fig. ^c). As indicated by the labels, in addition to the 
fundamental CDW peak (1), four more peaks lie along 
the CDW direction (2-5). Although the transforms of 
real space images resemble diffraction data, symmetry re- 
quirements intrinsic to diffraction data do not apply. As 
demonstrated by Thomson and co-workers, the Fourier 
transforms of STM images that exhibit true discommen- 
surations always have extra peaks in proximity to the fun- 
damental CDW peak . This arises from the fact that 
a discommensurate CDW can be expressed as the prod- 
uct of a uniform ly i ncommensurate CDW and a mod- 
ulation envelope [18j. The wavelengths of the envelope 



4 



i 1 1 1 1 1 1 ■ i 1 1 




■ (c) 



2 4 6 8 10 12 14 
bond index 

FIG. 4: Te-Te bond length deviation from the average value 
as refined crystallographically (a) and from the PDF (b). The 
deviation rbond — ro is defined such that n is the Te-Te bond 
length of the i th bond (bond-index i) in the unit cell and ro — 
3.1 A. (c) Schematic of the arrangements of "short" bonds 
within the unit-cell highlighting the formation of oligomers. 
Short bonds are defined as those lying within ± a of the center 
of the first Gaussian of the bimodal distribution in Fig. [3(b). 
indicated as a yellow band in the Figure. The Gray markers 
represent Te atoms. Red (blue) lines and markers are bonds 
lying in the top (bottom) row of the unit cell. The CDW is 
out of phase between these two rows of bonds. 



are given by the differences in the wave vectors of closely 
spaced peaks. The longest such wavelength in our im- 
ages is 38 A, corresponding to peaks 2-3 and 3-4 in the 
Fourier transform, indicating that a discommcnsuration 
separation of this length-scale exists. This is consistent 
with the refined PDF behavior which crosses over from 
the local to the average behavior for a refinement range 
of around 27 A(Fig. |3J), which would be expected to oc- 
cur at around, or a little above, the radius of the local 
domains. 

As well as the bond-length distributions, the local and 
average structure refinements allow us to study the pat- 
terns of atomic displacents due to the IC-CDW. The av- 
erage structure refinement |2lj results in an almost per- 
fectly sinusoidal pattern of bond-lengths, with the wave- 
length of the CDW (Fig. |2 a)), clearly identifying these 
as Peierls distortions. The local structural model was 
refined in the same unit cell and space-group, but re- 
sults in a much more square-wave like distribution, con- 
sistent with the distinct short and long Te-Te distances 
described above (Fig. Bib)). Fig. 0{c) shows the pat- 
tern of bonded Te-Te atoms that results when the short 
distances determined from the PDF data are plotted in 
the unit cell. In this way, the Peierls distortions due to 
the IC-CDW can be thought of as forming oligomers in 
the Te net. In this picture, the discommensurations oc- 
cur when the pattern of oligomers has defects. This is 



a common picture in the chemistry literature |2ll |30|, 
though we note that this picture is not supported by the 
crystallographic results shown in (Fig.^a)) and needed 
the application of a local structural method to show that 
it has a physical reality beyond its heuristic value. 

The refined parameters of the \ow-r max PDF re- 
finements yield quantitatively the atomic displacements 
within the commensurate domains. This is the first 
demonstration of the use of the PDF to obtain quan- 
titatively the atomic displacements (Peierls distortion) 
within the commensurate domains of a discommensu- 
rated IC-CDW. This opens the way to a quantitative 
first principles calculations and a better microscopic un- 
derstanding of the IC-CDW state. 
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